***************************************************
****----			Figure 4	  		   ----****
***************************************************

clear all
cap log close 
set more off

****---- Change directory

global data "/Users/Replication material/Data"

****---- Figure

use "$data/dataset_replication_regional.dta", clear


qui reg N_influenza_polmonitePC N_caduti_tot_frontePC  if anno==1918, r cluster(regione)
predict sign, res
qui reg N_influenza_polmonitePC N_caduti_malattiaPC  if anno==1918, r cluster(regione)
predict sign2, res
  
scatter sign N_caduti_malattiaPC  if anno==1918, xlabel(0.0001 0.0002 0.0003 0.0004, angle(45)) ylabel(, angle(45)) mlabel(regione) graphregion(color(white))  ytitle("Deaths - influenza (per-capita) [residual]") xtitle("Victims WWI - illness in home-town (per-capita)") legend(off) || lfit sign N_caduti_malattiaPC  if anno==1918

scatter sign2 N_caduti_tot_frontePC  if anno==1918, xlabel(, angle(45)) ylabel(, angle(45)) mlabel(regione) graphregion(color(white))  ytitle("Deaths - influenza (per-capita) [residual]") xtitle("Victims WWI - outside home-town (per-capita)") legend(off) || lfit sign2 N_caduti_tot_frontePC  if anno==1918


qui reg N_influenza_polmonitePC N_caduti_tot_frontePC  if anno==1917, r cluster(regione)
predict sign3, res
qui reg N_influenza_polmonitePC N_caduti_malattiaPC  if anno==1917, r cluster(regione)
predict sign4, res
  
scatter sign3 N_caduti_malattiaPC  if anno==1917, xlabel(0.00003 0.00005 0.00007 0.00009 0.00011 0.00013 0.00015, angle(45)) ylabel(, angle(45)) mlabel(regione) graphregion(color(white))  ytitle("Deaths - influenza (per-capita) [residual]") xtitle("Victims WWI - illness in home-town (per-capita)") legend(off) || lfit sign3 N_caduti_malattiaPC  if anno==1917

scatter sign4 N_caduti_tot_frontePC  if anno==1917, xlabel(, angle(45)) ylabel(, angle(45)) mlabel(regione) graphregion(color(white))  ytitle("Deaths - influenza (per-capita) [residual]") xtitle("Victims WWI - outside home-town (per-capita)") legend(off) || lfit sign4 N_caduti_tot_frontePC  if anno==1917
